Detecting Chiral Orbital Angular Momentum by Circular Dichroism ARPES 
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We show, by way of tight-binding and first-principles calculations, that a one-to-one correspon- 
dence between electron's crystal momentum k and non-zero orbital angular momentum (0AM) is 
a generic feature of surface bands. The 0AM forms a chiral structure in momentum space much 
as its spin counterpart in Rashba model does, as a consequence of the inherent inversion symmetry 
breaking at the surface but not of spin-orbit interaction. Circular dichroism (CD) angle-resolved 
photoemission (ARPES) experiment is an efficient way to detect this new order, and we derive 
formulas explicitly relating the CD- ARPES signal to the existence of 0AM in the band structure. 
The cases of degenerate p- and d-orbital bands are considered. 

PACS numbers: 



C/2 



X3 
O 

o 



00 



X 



Electron spins are quenched in ordinary crystalline 
solids in the sense that each crystal momentum k comes 
in degenerate pairs of spin-up and spin-dovifn electrons. 
Such degeneracy is lifted in an interesting manner for the 
so-called Rashba system[l], in which the given momen- 
tum state in the band has only one spin state associated 
with it. Examples of Rashba-split bands are many by 
now [2]. On symmetry grounds, as Rashba originally ar- 
gued, the inherent inversion symmetry breaking (ISB) at 
the surface termination allows an interaction term, the 
Rashba term, of the form iJ^ = Arz • (k x cr) involv- 
ing the coupling of the electron's spin operator cr/2 and 
its momentum k (We set h = 1). The chiral spin an- 
gular momentum (SAM) structure in momentum space 
follows as a direct consequence of the Rashba Hamilto- 
nian i?R[3, 4]. 

One must note, however, that other physical quan- 
tities having the character of angular momentum can 
take the place of cr in the Rashba Hamiltonian. For in- 
stance, in bands where the orbital angular momentum 
(0AM) remains unquenched, one can equally well con- 
sider a coupling of k to 0AM operator, L. In this paper, 
we expand this observation to show that chiral 0AM in 
one-to-one correspondence with the electron's linear mo- 
mentum is indeed a general consequence of ISB for the 
surface bands. Unlike the chiral SAM structure which 
typically occurs in materials with strong spin-orbit inter- 
action (SOI) [2], we show that chiral 0AM occurs even 
for materials with no SOI. 

We begin by providing a simple, tight-binding (TB) 
model Hamiltonian in two dimensions where the above- 
mentioned chiral 0AM structure emerges. A three-fold 
degenerate p-orbital system is considered, while the spin 
degrees of freedom is suppressed in order to emphasize 
the notion of chiral 0AM emerging from ISB alone with- 
out the requirement of SOI. The TB model is the triangu- 
lar lattice version of the square one originally considered 



by Petersen and Hedegard[5] in their study of Rashba 
phenomena in p-orbital bands. By suppressing spin de- 
grees of freedom, we do away with the notion of conven- 
tional Rashba spin splitting. The idea of ISB-induced 
chiral 0AM was not considered in Ref. 5 or in any other 
literature we are aware of [6] . 

We introduce two Slater-Koster parameters Vi and 
V2 for cr- and 7r-bonding amplitudes, respectively, and 
a third one, 7, representing the degree of ISB [5]. We 
will write N for the number of sites in the lattice, and 
|i. A) for the localized Wannier orbitals (A = Px^Py,Pz) at 
the atomic site r^. Then the tight-binding Hamiltonian 
for the triangular lattice (essentially the same result ob- 
tains for square lattice) near the P-point (k — 0) in the 



momentum-space basis |k. A) = 
comes 
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diagonalize H]^^, it is convenient to choose a new set of 
basis vectors 



fc^. The lattice constant is taken to be unity. To 
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k — |k|, e^'^'^ — {kx +iky)/k. The state |k, I) remains 
decoupled at energy Ei^ = 3^2 - 3^1 + 3(Vi - 3^2)1*2/8, 
while |k, II) and |k, HI) combine to form eigenstates 
which, to leading order of 7/ A, are 
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with energies E2,k ^ S{V2 - Vi) + 3((3Vi - V2)/8 ~ 
7V4A)k2 and Es^^ ^ &V2 + (S-f^/iA - ^¥2/2)]^^. The 
0AM operator is given by the sum L = {i/N)J2i'^i 
where each acts on the Wannier state |i, A) as the 
usual L — 1 angular momentum operator. The two bands 
obtained above carry nonzero, chiral 0AM as claimed 
(i+ = + iLy): 



(k,2|L+|k,2) = 2(fc, 



+iky) = — 



;k,3|L+|k,3). (4) 



The "helicities" of the chiral OAMs are opposite between 
the bands, with their magnitudes vanish linearly with |k| 
and the ISB parameter 7. In contrast the Rashba model 
displays perfect spin polarization irrespective of the wave 
vector or the strength of the Rashba parameter Ar. 

A second check on the existence of chiral 0AM is per- 
formed by employing the first-principles local-density ap- 
proximation (LDA) calculation for a Bi single layer form- 
ing a triangular lattice. The choice is inspired by Bi being 
a proto-typical p-orbital band material. An external elec- 
tric field of 3V/A perpendicular to the layer was imposed 
by hand to mimic the surface potential gradient without 
having the complication of dealing with the bulk states. 
We also performed calculations for the physically more 
realistic case of Bi bi-layer with perpendicular electric 
field [8], with results that are entirely in accord with the 
statements made below for the single-layer case regarding 
the emergence of chiral 0AM. To emphasize the relevance 
of ISB we again chose to investigate the spin-degenerate 
case by turning off SOI in the LDA calculation. The re- 
sulting electronic structure for spinless case consisting of 
three p-orbital-derived bands is shown in Fig. 1(a). As 
the external electric field is turned on, a level repulsion 
between the middle {E2 in Fig. 1(a)) and the bottom 
(£■3 in Fig. 1(a)) band occurs as indicated by circles in 
Fig. 1(a). These two bands exhibit the chiral 0AM pat- 
terns with the maximum 0AM vector \{L)\ w 0.96fi- as 
shown in Fig. 1(c) and 1(d), while the third one, shown 
in Fig. 1(b), carries much less 0AM around the F point. 
The 0AM chiralities of the two bands are opposite, in 
accordance with the previous TB analysis. An excellent 
fit of the LDA band structure near the F point was possi- 
ble with the TB parameters Vi = -0.725 eV, V2 = -0.11 
eV and 7 = 0.2623 eV (Fig. 1(a)). The 0AM magnitude 
is seen to decrease continuously upon approaching the 
F point in the LDA calculation (Fig. 1(c) and 1(d)) as 
predicted by the TB calculation, Eq. (4). 

Having established theoretically the existence of chi- 
ral 0AM in generic inversion-asymmetric bands by a 
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FIG. 1: 0AM and CD from first-principles and tight-binding 
calculations of Bi monolayer without SOI. (a) LDA band 
structure for Bi monolayer with SOI turned off. Perpendic- 
ular electric field of 3V/A was imposed externally. Three 
dashed curves represent the tight-binding energy dispersions 
around the F point, (b)-(d) OAM vectors (green arrows) and 
NCD signals (color backgrounds) for the three bands, i?i(b), 
E2{c), and -B3(d) over the whole Brillouin zone marked by a 
solid hexagon. Largest OAM has magnitude « Ih for bands 
E2 and E3. NCD is calculated with kF,z = 2.27A"^ (e)-(f) 
NCD calculated with kp.z = in Eq. (12) for E2{e) and £;3(f) 
bands. The opposite color assignments between (c) and (e) is 
a consequence of photon energy dependence of the scattering 
intensity. 



number of methods, we turn to the question of its de- 
tection. Spin- and angle-resolved photoemission spec- 
troscopy (SARPES) has served to identify the chiral spin 
structure of the surface bands in the past [3, 4]. A similar 
chiral structure for OAM as demonstrated here cannot, 
however, be detected by the same probe since chiral OAM 
exists even when SOI is very weak and spin degeneracy 
is nearly perfect. Given the potential prevalence of chi- 
ral OAM in ordinary, i.e. weak SOI, surface bands, it is 
essential to establish detection tools for this new degree 
of freedom in band materials and develop some theory of 
it. 

Circular dichroism (CD) refers to phenomena in which 
the physical response of a system to probing light de- 
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pends systematically on the light polarization being left- 
circularly-polarized (LCP) or right-circularly-polarized 
(RCP). In ARPES, CD manifests itself as different scat- 
tering intensities of the photo-electrons depending on the 
helicity of incident hght being RCP or LCP. The CD- 
ARPES signal can be quantified through the normalized 
CD (NCD) defined as 



Z?(k) = 



E.(/?CP(k)+jLCP(k))- 



(5) 



The sum over the final state spin a reflects the spin- 
integrated nature of the detection scheme. Spin index is 
restored in the following derivation of the NCD formula 
as it will apply equally well to strong-SOI materials with 
mixed spins. 

The initial state |I) is the Bloch state of momentum 
k constructed as |k, m) = iV~^/^ e*'''"'' |i, m). The 
Wannier state \i, m) at site r, is given by 



vector G. Making use of the fact that /(|r|) depends only 
on the radial distance, one can re-write the second line 
of Eq. (8) as 

-A^m^,„dc,d^ J e-'''^--/(|r|)d^r = -A • Vk^ff<,(kf ),(9) 



where we have introduced 



-I' 



fi\r\)d\ 



(10) 



and = {mx,cr,'my,cr,'mz,cr)- One finds the average 
0AM (L) for the p-orbital Bloch states |k, m) expressed 
simply in terms of mo-: 

(k, m|L|k, m) = (L) = i ^ m<,(k) x m;(k). (11) 



,m) = ^mA,a(k)|i,A,a), 



(6) 



as a linear combination of constituent atomic orbitals 
labeled by A, and spin a, with k-dependent coefficients 
mA,<T(k). Usual plane wave forms are assumed for the 
final state: Vf ~ e^"". 

The transition amplitude into the final state of spin a 
is evaluated as 



(F,a|p-A|/)^(F,a|r-A|/) 

~^m;,,,(k)(F|r.A|i,A)e''^-'-v 



(7) 



Given the localized nature of the Wannier state, it is 
useful to re-write r • A = (r — r^) • A -I- • A. One 
immediately finds that {F\ri • A|/) = (r^ • A)(F|/) is zero 
from the presumed orthogonality of the initial and the 
final states. To proceed further, we treat the cases of 
p-orbital and c?-orbital bands separately as they require 
somewhat different strategies for evaluation of the NCD 
formula. For p-orbitals the Wannier states are assigned 
the hydrogenic wave function (r|z, A) ^ (r — ri)\f(\r — 
Til) (A = x,y,z) and with these we find the transition 
amplitude becomes 



X^mA,.(k) J e"*'^-'-[r.A]xA/(|r|)j 



-ikj?-r I 



r-A]xxf{\r\). 



(8) 



The sum e'^^-^'^l-"^ = (5(k|,-k+G) simply yields in- 
plane momentum conservation up to the reciprocal lattice 



Of great importance is the fact that nonzero 0AM is 
possible even with spin degeneracy, rria = m, where one 
would obtain (L) = 2im x m*. We can now proceed to 
show that NCD is fundamentally related to (L). Equa- 
tions (8) and (10) allow the NCD formula (5) to be re-cast 
in the compact, suggestive form: 



D{\.) = 



(A X A*) • V,g, X Vgl 



[(A • Vg.) (A* ■Vg*) + (A* • (A- V.9* )] 



^2) 



Explicit reference to kp in the derivatives has been 
dropped from above. The vector potentials are A = 
{ei+i£2)/y/2 for RCP and A* for LCP, with A x A* = 
—isi X £2 = — jfcph giving the photon direction. The re- 
maining task is to evaluate the quantity J2a- ^ 
which governs the CD response of the given initial state. 
The denominator, by definition, is always positive and 
plays a minor role in characterizing the 0AM structure. 
For p-orbitals, inserting g„ = ■ V/ gives 



^ X Vgl = le'^^'^im, x m;)„Va^/ x Vd^f 



= -_e'-I^^L^)^Q^f X Vd^f, (13) 

where Eq. (11) has been adopted in the second line. 

Clearly, this is proportional to the components of 0AM 
in the initial state. The remaining "form factor" 'Vdjsf x 
V9-y/ depends onkp, which in turn depends on the in- 
coming photon energy. Both (L) and the form factors 
can be obtained by faithful LDA calculations of the wave 
functions. Analytically, choosing /(r) e"""/", for in- 
stance, yields f{kp) = 1/ {l+kpO,^)^ and the form factors 
can be worked out. 
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Next we turn to the case of degenerate rf-orbital bands. 
Now g„ = (rn„ ■ D)/ involves the inner product between 
the five-dimensional coefficient vector m„ and the cor- 
responding differential operators matching the given 
orbital basis. Examples are D^y = d^dy, D^^2_^2 = 
{2dl-dl-dl) /2\/3, etc. Although the dichroism formula 
(12) still applies, we arc no longer able to transform its 
numerator to a simple form like Eq. (13). Recall that in 
realistic ARPES experiment the incident photons carry 
energies of several tens of cV, delivering however at most 
an eV of energy to the occupied electrons [7]. After sub- 
tracting what amounts to the surface potential energy 
barrier, there is still a lot of energy imparted to the fi- 
nal photo-electron, whose energy is typically in excess of 
10 eV. Due to in-plane momentum conservation, the in- 
planc component of photo-electron momentum can carry 
only a small fraction of this energy, which means most 
of the kinetic energy is contained in the ^;-component, 
(kp.z)'^ /2m. As a result, the typical situation in ARPES 
experiment is the one in which kF,z dominates over the 
planar components, and kF,za is rather larger than unity, 



a being the radius of the orbital wave function. 

In computing g„ = {rna- ■ D)/, therefore, the orbitals 
containing at least one power of z will be dominant over 
those that contain none, due to extra powers of kpU 
produced by the differentiation. The reasoning reduces 
the number of relevant orbitals from five to three, i.e. 
dzx, dyz, and (13^2 _r2. This is irrespective of particular 
crystal field symmetries of the d-orbitals and is rather 
dictated by the experimental conditions of ARPES. We 
also introduce the "partial 0AM" obtained by assuming 
only the three, labeled 1 = zx, 2 = yz, and 3 = 3z^ — r^, 
out of the five d-orbitals contribute to the wave function: 
{L^y = VSiJ2a{"^2,aml^ - ml„m3,ay {LyY = 



The prime is a re- 



minder that contributions average 0AM from d^y and 

d^2_y2 are missing. Following the same calculation 
procedure as in the p-orbital case, we find the CD 
formula for d-orbitals valid for kF,zO' '> 1, 



fcph • {2{Lxyx + 2{Ly)'y - {Lz)'z) 



Ea[^xmi,<, + Aym2,a - 2y/2,A,rm,^][{A^)*ml^„ + {AyYml^^ - 2s/^{A,)*ml„] + {A o A*) " 



(14) 



The 0AM shown in the numerator refers to those in partial 0AM. Again, non-trivial CD- ARPES signal is predicated 
on the existence of 0AM. In the same limit, fcir.zO ^ 1, the p-orbital formula simplifies to 



5fcph • (L)-6(fcph • kF){kF ■ (L)) 



Y.MkF ■ A){kF ■ m„) - {m„ ■ A)][6(/cf • A*){kF ■ m%) - (m* • A*)] + (A ^ A*) 

I 



(15) 



The case of <2g-bands involving the three orbitals, yz = 1, 
zx = 2, and xy = 3 can be considered as well. The factor 
E(7 ^5<T X V^; in Eq. (12) becomes 



(16) 



where D = {Dy^, D^x, D^y). As in the p-orbital case, we 
obtain a simple relation for the 0AA/[ average of the given 
Bloch state in the i2g-band: (L) = —iY^^fria x m*. 
(Note the opposite sign compared to the p-orbital result 
in Eq. (11).) Again, non-trivial 0AM is responsible for 
NCD for t2g-bands. The actual task of evaluating the 
form factors leads to cumbersome expressions. For prac- 
tical applications, it is better to evaluate them numeri- 
cally using LDA-obtained wave functions. 

The 0AM of the Bi bands shown in Fig. 1(a) is used to 
obtain D{k) according to Eqs. (12) and (13). Exponen- 
tial function /(r) ~ e"*"/" with the radius a = 1.6Awas 



used to evaluate the form factors. Figure 1(b)- (d) shows 
D{k) overlaid with local 0AM when = 2.27A. The 
next set of figures in (e) and (f) is obtained with kp- = 
OAfor E2 and £^3 bands. Incident photon direction is cho- 
sen fcph = (cos 60°, 0, — sin 60°) in both sets. The photon 
energy dependence of the NCD is obvious by comparing 
Fig. 1(c) with (e). LDA calculation for light-element 
d-orbital bands are under way and will be reported else- 
where, along with some experimental results [9]. We com- 
ment that a potential modification of the NCD formula 
may be necessary for heavy elements due to significant 
spin-orbit interaction in the Hamiltonian and the conse- 
quent modification of the current (momentum) operator 
coupling to A[10]. Our present results are, however, free 
from such complications as we exclusively focus on spin- 
degenerate band structures. Non-trivial 0AM and NCD 
features are only due to the degenerate multiple-orbital 
nature of the eigenstates. 
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